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1. Introduction 

"Natura abhorret vaccum" was an antique saying. Modern science has strangely supported 
this intuition when discovering that vacuum is very far from being empty. There are quantum 
fluctuations, virtual e^e^ and qq pairs and vacuum expectation values (v.e.v). Our target here is 
the v.e.v of = AaA"^ in Landau gauge. 

Is it legitimate to speak about a gauge dependent property of a gauge invariant vacuum ? 
Yes, of course, once we have understood that there are fields in the vacuum. These fields may 
be scrutinized in different gauges, the appearance will differ, the physics we learn will vary, and 
nevertheless the vacuum is gauge invariant. It is absolutely similar to the well known fact that a 
physical process, say some reaction, is decomposed differently into Feynman diagrams according 
to the chosen gauge, its appearance will depend on the gauge although the process is gauge invariant 
and the cross section will not depend on the gauge. 

A^ is certainly not a gauge invariant quantity, but in Landau gauge it is invariant for infinitesi- 
mal gauge tranformations, and consequently, for BRST transformations. 

Why do we care about (A^) ? It is useful 

• in order to understand the infrared properties of QCD. It is enough to quote an incomplete 
list of studies of QCD in the infrared which use (A^) [1]-[15], many of the authors being in 
the audience; 

• in order to identify correctly non-perturbative corrections to renormalisation constants; 

• in order to test and calibrate on a rather extensively computable case Operator Product Ex- 
pansion : we find larger non-perturbative corrections than usually expected. 

2. How to define (A^) ? 

A naive estimate of (A^) produces an ultraviolet divergent quantity oc {a being the lat- 
tice spacing). Indeed we are interested in infrared modes and we must define a scale jU which, 
grossly speaking, separates the high and low modes. A way to perform rigorously the distinction 
between UV and IR modes is to use Operator Product Expansion (OPE). One must then define 
a renormalisation scheme and a renormalisation scale. 

3. Wilson Operator expansion 

A momentum dependent quantity Q{p^) with vacuum quantum numbers can be expanded in 
inverse powers of /i^ 

!2(/) = !2pert(p^At') +ce,^„„(/,M')(A'(At')) + •••• (3.1) 

where 2pert(;?^,M^) and C^iison(/'^)M^) ~ series in o;(/x^) computable in perturbative 

QCD. Q can be the strong coupling constant, the quark field renormalisation constant, other renor- 
malisation constants, etc. The coefficients C^jisQ„ are often called Wilson coefficients. 

OPE has been extensively used in phenomenology since the pioneering work by SVZ [16]^. 

'And we have the pleasure to have "Z" in the audience. 
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3.1 Criteria to check that we really measure 

• The use of OPE is criticized under the argument that it is difficuh to distinguish higher 
orders in the perturbative series (which are logarithmically suppressed oc l/log(/7^)) from 
OPE contributions which are power suppressed (o^ 1/ p^) and that perturbative renormalons 
mimic a condensate, being also oc 1/ p^. This issue has been carefully studied by Martinelli 
and Sachrajda (MS) in [17]. They show, using a specific model, that when subtracting with 
proper factors two quantities, say in our case Zq and a^, renormalon ambiguities cancel. One 
must still check that the renormalon-free neglected terms in the perturbative expansion are 
not larger than the condensate. The criterium proposed by MS is that we must check that 
the last perturbative contribution considered (say the order or depending on what is 
available) is significantly smaller than the power connection. This criterium has been checked 
in [14] concluding that the MS rule was satisfied on the full momentum range which is used. 

• is not a lattice artefact. It is a quantity defined in the continuum limit. Therefore 
we must check that what we interpret as {A^{pL^)) does not depend significanlty on the lattice 
spacing. 

• We must also check that (A^(/x^)) derived from different quantities are consistent provided 
we use the same renormalisation scheme and scale. 

• {A^{n^)) depends on the vacuum and there is no reason for the values of {A^{jj})) extracted 
from Nf = and Nf = 2 lattice data to agree. 

3.2 About Wilson coefficients 

All Wilson coefficients C^iison(P^iM^) of the type 



C^nson(/,M') =^^treeg^(M^) '-^^ where 



p2 



1 9 

t/2tree=— for 2 = Z,,, d^tvee = —foi: Q = aT. (3.2) 
12 32 



From eq. (3.1) and eq. (3.2) we see that there is always the same factor and 
we will therefore give the fitted values of the condensate as {g^{pL^)A^{ll^))y^^. Indeed, in the 
following, we choose to renormalise the Wilson coefficient such that the local operator A^(/i^) is 
in the MS scheme, whichever prescription we take for the perturbative part of Q. 
In practice one can show that the best and most general fitting formula is: 

e(/) = Gpen(/,M') + ^Q^°jpl'^l^^ {A\n^))u)j (3-3) 
At leading logarithm for the non-perturbative correction, 

''^"^""^''''^'^ (A^(M^))MS = ^^tree {g\^')A\n'))^ 



44 



where e=— — (3.4) 
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Notice that e is small. Therefore the corrective factor in eq. (3.4) is almost scale invariant. The 
exponent is the same for all quantities. In eq. (3.4) the only term which depends on the measured 
quantity is cf^tree. 

3.3 How to compute (A^(/x^))f^s from LQCD ? 

• Compute 2iatt as measured from LQCD. 

• Correct for hypercubic lattice artefacts i.e. artefacts related to the hypercubic geometry. 

• Fit according to 

where Caipi'^^P^ is a simple model for the non-hypercubic remaining lattice artefacts, which 
turns out to give a good result, as we shall illusti^ate below. 



Let us now apply this strategy. 



4. The strong coupling constant 

There are many ways to define the strong coupling constant. We will use [18, 22] what we call 
the "Taylor coupling constant" which, thanks to Taylor's theorem [19], is only dependent on the 
gluon and ghost propagators. 

We will use configurations with Wilson twisted quarks (Nf=2) from the ETM collaboration 
and compare them to quenched configurations (Nf=0). 

4.1 Some definitions 

The gluon propagator (G^^^) and the ghost propagator (F^^^) are defined as follows, G{p^,K) 
and F{p^,A) being named the "dressing functions": 

{f^^^)'\p\K) = -sJ-^; (4.1) 

A being some regularisation parameter: A = (/3) if, for instance, we specialise to lattice regu- 
larisation. The renormalised dressing functions. Or and Fr are defined through : 



Gr{p\^^) = limA^^Z3-njU',A) G(p2,A) 

FRip\^^) =\imf,^^Z-\pL\K)F{p\K) , (4.2) 
4 



(A^) from LQCD in Landau gauge 



O. Pene 



with renormalisation condition 

G«(Ai',At')=fR(At',M') = l • (4.3) 

Now, we will consider the ghost-gluon vertex which could be non-perturbatively obtained through 
a three-point Green function, defined by two ghost and one gluon fields, with amputated legs after 
dividing by two ghost and one gluon propagators. This vertex can be written quite generally as: 



rf'i-g,k-q-k)=igof''{qyH,{q,k) + {q-k)yH2{q,k)) , (4.4) 
The vertex renormalisation constant is defined as 

{Hf{q,k)+Hi{q,k))L , = limZl(A^^A) {Hi{q,k;A)+H2{q,k;A))\ 2^ 2 = 1, (4.5) 
The renormalised coupling constant is defined as 

gR(At') = lim Z3(A^^A)z]/'(Ai^A)^o(A') (H,{q,t,A)+H2{q,t,A) 

= hm go{A') ' I ' . (4.6) 

Now we choose a special kinematics where the incoming ghost momentum vanishes. Taylor's 
theorem states that//i (<7,0;A) +H2{q,0;A) is equal to 1 in full QCD for any value of q. Therefore, 
the renormalisation condition eq. (4.5) implies Zi(jU^) = 1 and then 

ar(At') = = lim ^^Gi^\A^)F^{^\A^) ; (4.7) 

471 A^oo 471 

which only depends on the propagators. 

4.2 Perturbative running of aj 

The four-loops expression for the coupling constant in the Taylor scheme as a function of Aj 
(Aqcd in this scheme) is given by [20, 21, 22]: 

2 

where ; = In ^ and coefficients are 

/3o = 11 -3A^/, i3i=jSi = 102-yA^/ 
^2 = 3040.48 - 625.387 iV/ + 19.3833 A^^ 

^3 = 100541 - 24423.3 Af/ + 1625.4^) - 27.493^) , (4.9) 
At is converted into A^ by 

ci 507 -40A^/ 

^ = ~Wo = 792 -48^/ (4 
Ar 
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4.3 The non-perturbative contribution to ocj- 



It is easy to see that the dominant non-perturbative contribution to ar(/x^) is the condensate 
(A^). To take it into account in our fits we will need the Wilson coefficient C^fj^^j^. This is possible 
up to three loops thanks to [8] where the (A^) correction to the ghost and gluon propagators have 
been computed. We will only write here the result at leading logarithm [22]. 



An^ 



-9/(44-8^/73) 



1 + 



^QCD 



s\{pil){A')R. 



\ 



In 



\ \ ^QCD , 

where /Xq is our reference renormalisation scale which we take to be 10 GeV. 



(4.11) 



4.4 Results of the lattice simulation 




We perform the lattice calculation of 
arilJ.^) from eq. (4.7), using Nf = 2 
with twisted dynamical quarks from the 
ETMC collaboration. We eliminate the 
hypercubic lattice artefacts. The results 
depend on the dynamical quark mass. 
For three values of the lattice spacing 
(labelled by j8) we extrapolate to zero 
mass leading to the three sets of data in 
the upper-left plot of the figure aside. 



We invert eq. (4.8) to extract Aj- as a function of /i^ and then Ay^ from eq. (4.10). Since Aj- and 
A^ are constants independent on the scale, the result should be a nice "plateau" if we were in 
the full perturbative regime. This is obviously not the case as shown in black for the three lattice 
spacings in the upper-right and two down plots. We then use eq. (4.11), fitting g^(jU^)(A^)^ ^2 
to have A^ as a nice plateau. The result is shown in blue. It looks quite convincing. From the 
comparison of the black and blue curves it is evident that, although the fit has been performed 
in the rather high-energy range 2.6-6.4, a sizeable non perturbative correction is needed. 
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0) 




=0) 




0) 


a^; a A^O 


134(7), a^g^A 


--■ aj;aA=0 


134(7), a"gi:A 



a(P=3.9) p 



In the plot on the left we have merged the three 
curves of arin^) by fitting the lattice spacing 
ratios to have the best matching. The compati- 
bility of the three curves is quite impressive and 
the needed lattice spacing ratios agree very well 
with what is extracted from the interquark poten- 
tial [22]. The fit has been performed both us- 
ing the leading logarithm [22] (LL) and the three 
loop [8] (^(a"^)) formula for the Wilson coeffi- 
cient. The result is: 



Ams = 330±32+°3 MeV 



(4.12) 

g^i^^){A^)R,^ = 2.7± 1.0±0.7GeV2(^(a4)) 



g\lx^){A^)R,^ = 4.4 ± 1.5 ±0.7 GeV2(LL) 

5. The quark field renormalisation constant Zq 

• We compute the Fourier transform of the lattice quark propagator for momentum p: it is a 
12 X 12 matrix S{p). We define the renormalised quark field 



^R=Zg^bare whcnCC Zg{n =p) 



12p2 



rTr 



(5.1) 



We suppress the hypercubic artefacts which in this case ai^e particularly large. 
We perform the fit 



^qlatt 



b^) =Zpe«(/,M2) + +^a2p2aV (5-2) 



4 



.9 -< ratio i — i — 
.8 -< ratio < .9 i------- 

.7 ^< ratio< .8 
.6 -< ratio < .7 □ 
.5 -< ratio < .6 
.4 -< ratio < .5 
.3 -< ratio < .4 ;- 
.2-<ratio<.3 - - 
, hyp corrected OWF h---'--- 



H. ■ 



The plot on the left shows the re- 
sult for Zqiatt as a function of p^ 
stemming directly from the lattice 
calculation. It is very fai" from 
the smooth dependence expected 
in the continuum. One sees a 
"half-fishbone" structure which is 
a dramatic expression of hypercu- 
bic artefacts. Let us define 



J4] 



ratio ; 



(p2)2 
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In the above plot the color code corresponds to the value of the parameter "ratio" which is bounded 
0.25 < ratio < 1. It is visible that the hypercubic artefacts increase with this parameter. We reduce 
drastically these artefacts by an extrapolation down to ratio = [12] : 

n[4] 

+ Ca4p4a^ (5.3) 




The merged plot results with j3 = 4.05 and j3 = 4.2 rescaled to the j3 ~ 3.9. The l.h.s shows the data 
corrected for all lattice artefacts. The r.h.s shows the same data furthermore corrected by the perturbative 
running factor up to 10 GeV. The horizontal axix is in GeV ^. The black line on the l.h.s corresponds to 
the global fit with perturbative running and the three-loops (Chetyrkin-Maier [8]) Wilson coefficient for the 
1 / term. The black Une on the r.h.s corresponds only to the 1 / times the three loops wilson coefficient 
added to Zf ((10 Ge V)^ , 6/3 .9) = 0.726 

Once this non-perturbative hypercubic correction is performed, we merge the results for three 
lattice spacings. The merged result is shown in the figure above l.h.s. In the r.h.s these data are 
corrected by the perturbative running factor up to 10 GeV. If non-perturbative coiTcctions were 
absent the curve should be flat. This is obviously not the case. The results are shown in table 5. 





a^fm^ 










3.9 


0.00689 


0.726(5) 


0.0201(13) 


3.20(38) 


2.62(31) 


4.05 


0.00456 


0.742(5) 


0.0200(15) 


3.09(65) 


2.57(54) 


4.2 


0.00303 


0.760(3) 


0.0194(8) 


3.23(55) 


2.74(47) 


average 






0.0201(3) 


3.18(28) 


2.64(23) 



Table 1: Results for TF^^^ (lOGeV) and c^api eq. (5.2) and the estimated g^(A^) v.e.v from the tree level 
1/p^ term and from the Chetyrkin-Maier [8] (CM) Wilson coefficient. 

6. conclusion 

It is striking that the fitted value of g^(A^) are independent of the lattice spacings. It is even 
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more striking that the values extracted from the quark renormaUsation constant and the couphng 
constant are perfectly compatible. 

This is very encouraging in the sense that all our criteria to identify really a g^{A^) condensate 
are fulfilled. 



Nf 


order ^2 ^^2 


2. 


aj 


3 gluons 





LL 


9.4(3) 


5.2(1.1) 


10(3) 







9.0(3) 


3.7(8) 




2 


LL 


2.7(4) 


4.4(1.6) 




2 


0{a^) 


2.55(36) 


2.7(1.0) 





Table 2: Comparison of estimates of ^^(A^) from different quantities alNf — Q and Nf — 2. All are taken at 
the scale jx = \0 GeV. LL means leading logarithm for the Wilson coefficient. 0{a'^) refers to Chetyrkin- 
Maier's [8] computation. 

In table 2 we perform a global comparison of the estimates of g^(/x^)(A^)^ both in the quenched 
and unquenched case. The agreement between g^(/i^)(A^)^ estimated from different observables 
is not very good in the quenched case. This does not induce in our mind any doubt about the ex- 
istence of a g^(/i^)(A^)^ {pL = 10 GeV) condensate in Landau gauge for the following reason: we 
have performed a large number of fits with different inputs, and we have never found g'^{lj}){A^) ^ 
compatible with zero at better than four sigmas. On the other hand the fitting procedure is certainly 
delicate due to many correlations between g^(jU^)(A^)^ and lattice artefacts. Our final conclusion 
is that the accurate estimate of its value needs some improvements in our fitting method due 
to several correlations difficult to disentangle. There are however strong evidences in favor of 
the existence of a positive condensate g^(/i^)(A^)ju in the range 2-10 GeV^ in the MS scheme 
at = 10 GeV. 
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